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Abstract 

We are interested in prey-predator communities where the predator population evolves 
much faster than the prey’s (e.g. insect-tree communities). We introduce a piecewise 
deterministic model for these prey-predator communities that arises as a limit of a micro¬ 
scopic model when the number of predators goes to infinity. We prove that the process 
has a unique invariant probability measure and that it is exponentially ergodic. Further 
on, we rescale the predator dynamics in order to model predators of smaller size. This 
slow-fast system converges to a community process in which the prey dynamics is averaged 
on the predator equilibria. This averaged process admits an invariant probability measure 
which can be computed explicitly. We use numerical simulations to study the convergence 
of the invariant probability measures of the rescaled processes. 
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1 Introduction 

Prey-predator communities represent elementary blocks of complex ecological communities 
and their dynamics has been widely studied. The coupled dynamics of the prey and preda¬ 
tor populations is often described as a coupled system of differential equations. The most 
famous of them was introduced by Lotka |34| and Volterra | l5l in the 1920’s. There exist also 
stochastic models for these prey-predator communities as coupled birth and death processes 
(see Costa and al. m or as stochastic perturbations of deterministic systems (e.g. Rudnicki 
and Pichor [43]). 

In this paper we are interested in prey-predator communities in which the predator dynam¬ 
ics is much faster than the prey one. Such communities are common in the wild, especially if 
we consider the interaction between trees and insects (see Robinson and al. |42j for the study 
of Aspen canopy and its arthropod community or Ludwig and al. |35| for the interaction 
between spruce budworm and the forest). In these communities, the number of predators 
is much larger than the prey number and the predator mass is smaller than the prey one. 
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As a consequence, the reproduction and death events will be more frequent in the predator 
population than in the prey population. Such slow-fast scales have been studied in some 
Lotka-Volterra systems, mainly in the case of periodic solutions (e.g. HU)- In the following, 
we study two successive scaling (large predator population and small predator mass) of a 
microscopic model of the prey-predator community and introduce new stochastic processes 
for slow-fast prey-predator systems which corresponds to these scalings limits. In particular 
one interest of our work is that these processes never assume that the prey population size is 
infinite, as it is the case for models based on ordinary differential equations. 

We introduce a hybrid model for the demographic dynamics of the community where 
the prey population evolves according to a birth and death process while the dynamics of 
predators is driven by a differential equation. The community has a deterministic dynamics 
between the jumps of the prey population. This piecewise deterministic process arises as 
limit of a prey-predator birth and death process, when the number of predators tends to 
infinity while the prey number remains finite. Such piecewise deterministic Markov processes 
(PDMP in short) where introduced by Davis in 1984 [22][23]- They are used to model different 
biological phenomena. As an example, the dynamics of chemostats has been described as a 
piecewise deterministic model [201 DU ES E21- Chemostats, in which bacteria evolve in an 
environment with controlled resources, correspond to the opposite setting where the prey 
population (the resources) evolves faster than their predators (the bacteria). Other examples 
can be found in neuroscience, to model the dynamics of electric potentials in neurons (see 
Austin 0) or in molecular biology, where piecewise deterministic processes appear as various 
limits of individual based models of gene regulatory networks when the different interactions 
happen on different time scales (see Crudu and al. uni). 

In this paper we study the long time behavior of the prey-predator community process. A 
vast literature concerns the long time behavior of continuous time Markov processes. In the 
setting of piecewise deterministic processes, general results have been obtained by Dufour and 
Costa on the relationships between the stationary behavior of the process and a sampled chain 
(see [251 18]). We focus on the theory of Harris-recurrent processes (see Meyn and Tweedie 
[38 1 139] and references therein) that relies on Foster-Lyapunov inequalities. These inequalities 
satisfied by the infinitesimal generator of the process, ensure that the populations do not 
explode in some sense. Combined with irreducibility properties, they ensure the existence of 
a unique invariant probability measure and that the semi-group of the process converges at 
exponential rate to this measure. The irreducibility of the process is non trivial because the 
randomness only derives from the jumps of the slow component. Our proof relies on a fine 
analysis of the trajectories of the process. 

Further on, we rescale the predator dynamics by dividing the coefficients of the predator 
differential equation by a small parameter e. This scaling derives from the metabolic theory 
and illustrates the fact that the predator mass goes to 0 while the prey mass remains constant. 
The metabolic theory links the mass of individuals with their metabolic rates. Numerous 
experimental studies display relationships between the individual mass and the birth and 
death rates or the community carrying capacity (see Brown and al. uni, Damuth |2lj). Here, 
we simplify these relationships by assuming that the predator metabolic rates increase as the 
invert of their mass. This slow-fast system converges as e goes to 0 to an averaged process. 
In the averaged community, the predator population will always be at an equilibrium that 
depends on the prey number. Therefore the prey population evolves as a birth and death 
process where the predator impact is constant between jumps. In this case, computations 
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concerning the stationary behavior of the averaged process are easier because the community 
is fully described by the discrete dynamics of the prey population. 

This paper is organized as follows. In section [2] we present the piecewise deterministic 
process and the main results of this article. We give the first properties of the piecewise 
deterministic prey-predator community process and explain how it derives from an individual 
based prey-predator community process. In section [3] we study the ergodic properties of the 
prey-predator community process. These properties derive from a Foster-Lyapunov inequality 
and the irreducibility of the continuous time process and of discrete time samples. In section 
[4] we rescale the dynamics of predators and prove the convergence of the slow-fast prey- 
predator community to the averaged process. We prove that this averaged community admits 
an invariant distribution and study the convergence of the sequence of invariant measures of 
the slow-fast process as e —> 0 with numerical simulations. Finally we discuss in section [5] our 
results in view of biological and ecological applications. 

2 Model and main results 

2.1 The piecewise deterministic model 

We consider a community of prey individuals and predators in which the predator dynamics 
is faster than the prey dynamics. The community is described at any time by a vector 
Zt = (Nt, Ht) where Nt £ N is the number of living prey individuals at time t and Ht £ M+ 
is the density of predators. 

We assume that the prey population evolves according to a birth and death process. The 
individual birth rate is denoted by b > 0, the individual death rate by d > 0. The logistic 
competition among the prey population is represented by a parameter c > 0. The predation 
intensity exerted at time t on each prey individual is BHf. 

The predators density follows a deterministic differential equation whose parameters de¬ 
pend on the prey population. The individual birth rate at time t is rBNt. It is proportional 
to the amount of prey consumed by the predator. The parameter r £ (0,1) represents the 
conversion efficiency of prey biomass into predator biomass. The predator individual death 
rate D + CHt includes logistic competition among predators (D > 0, C > 0). 

The community dynamics is given by the differential equation 

j t H t = H t {rBN t -D-CH t ), (1) 

coupled with the jump mechanism 

P {N t+S = j\N t = n, H t ) = bns + o(s) if j = n + 1, n > 1, 

= n(d + cn + BHt)s + o(s ) if j = n — 1, n > 2, 

= 1 — (b + d + cn + BHt)ns + o(s ) if n = j. n > 2, (2) 

= 1 — bs + o(s) if n = j = 1, 

= 0 otherwise. 

Between the jumps of the prey population process N, the dynamics of the predator density 
H is deterministic. If the process Zt is at a point ( n, h) £ N* x M + after a jump, then the 


3 


predator dynamics is governed by the flow </ n associated with equation 0. More precisely, 
6 n satisfies: 


d_ 

dt 


4>n(h,t) 


4> n (h, t)(rBn - D - C(j) n (h, t)), fi n (h, 0) 


h. 


(3) 


Then for all i > 0, 


4> n {h,t) = 






(4) 


For h >0, the solution (f> n (h,t) remains positive for all t > 0 and converges as t —> oo toward 
an equilibrium h* given by 


h 


* 

n 


rBn — D 
C 


VO, 


(5) 


where a V 6 stands for the maximum of a and b. For sake of simplicity we introduce the global 
flow on N* X M_|_ h),t) = (to, <j) n (h, t )) for (to, h ) 6N*x M + and t > 0. 


In the following, the state space of the prey-predator process is denoted by E = N* x M + , 
we also define the subset E' = N* x (/i|,+oo). A generic point z G E is a vector ( n,h ) 
with n E N* and h E M+. The process Zt = (Nt, Ht)t>o belongs to the class of Piecewise 
Deterministic Markov Processes introduced by Davis (see |23|). It is a if-valued Markov 
process whose infinitesimal generator 


Af(n , h ) = h(rBn — D — Ch)d2f(n, h) + ^/(n + 1, h) — /(to, hbn 
+ ~ 1, h) ~ f(n, h)^jn(d + cn + Bh)l n > 2 . 


( 6 ) 


is well defined for functions / : N* x M + —> M bounded measurable, continuously differentiable 
with respect to their second variable with bounded derivative. The domain of the extended 
generator ([6]) has been characterized by Davis (Theorem 26.14 in {23]). 

We denote by Pt the transition semi-group and by (or P z ) the law of the process with 

initial condition z = (n,h) E E. 


Remark 2.1. In this model, we assume that the prey population cannot become extinct since 
the death rate is 0 when there is only one prey individual left. When this assumption is not 
satisfied, the prey population process can be dominated by a population process without predator 
which evolves as a logistic birth and death process. It is thus absorbed in 0 in finite time. 
The non-extinction assumption for the prey population is biologically relevant for trees-insects 
communities for example where the tree population rarely disappears thanks to a migration 
of trees (e.g. seeds driven by the wind). Here we chose to replace the migration probability 
of new prey individuals by the non extinction of the prey population. This choice allows to 
describe the prey-predator dynamics only with individual metabolic parameters such as birth 
and death rates. However, it is possible to include explicitly a migration at rate m > 0 in the 
prey population. We therefore define an alternative process (N^ m \H^ m ^) whose infinitesimal 
generator is given by 


A^f(n, h) = h(rBn — D — Ch)d 2 f(n, h ) + ^/(n + 1, h) — /(to, h)^j (bn + rro) 
+ ^/(to — 1 ,h) — /(ra, h) S jn(d + cn + Bh). 


(7) 
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The process H^) is a piecewise deterministic Markov process taking its value in 

_E' m ) = E U {0} x M + . In the sequel, we will mention when our results hold for this al¬ 
ternative model including migration and explain the modification induced by the migration. 


2.2 Main results 


Our main questions on the prey-predator process are twofold. First we are interested in the 
long time behavior of the process Z. 

Theorem 2.2. The community process (Zt)t> o is exponentially ergodic. It converges toward 
its unique invariant probability measure it at an exponential rate. There exist 0 < p < 1 and 
0 < R < oo such that, for all z = (n, h) E E', 

|| Pt(z, .) - 7r|| tv < Rp t (n + hr _1 ), 


Section [3] is devoted to the proof of the exponential ergodicity of Z. Our theorem relies 

In 


on the theory of Harris recurrent processes, whose main results are recalled in Section 3.1 


the setting of piecewise deterministic processes these results have been used to derive ergodic 
properties of different processes (e.g. an additive increase multiplicative decrease process |30j. 
a stress release process [33] . or a wealth-employment process m)- There exists also general 
results proven by |8| in the specific case where the deterministic dynamics admits a compact 
positive invariant set. The case of the prey-predator process Z = (N, H ) is more complex since 
the process is in dimension 2 and neither the jump rates nor the deterministic trajectories are 
bounded. 


Our second main result concerns a scaling limit of the process Z corresponding to the 
biological assumption that the predator mass is small. We introduce a sequence of processes 
Z £ for the rescaled parameters r £ = r/e, D £ = D/e and C £ = C/e. As e tends to 0, the 
dynamics of the fast component H £ is accelerated between the jumps of the slow component 
N £ . Therefore we expect that in the slow-fast limit, the prey population only depends on the 
equilibrium /i* of the predator density. The following theorem is a simplified version of our 


convergence result stated in Theorem 4.2 


Theorem 2.3. Fix T > 0 and assume rB — D > 0. We suppose that the sequence of initial 
conditions (Zq)o< £ <i converges to Zq in law and moreover that 


sup E((A r g) 4 ) < oo, sup E((Hq) 4 ) < oo. 
0<£<1 0<£<1 


Then the sequence N £ converges in law toward N in D([0,T],N) as e —> 0. 

The process N is a pure jump process on N* whose infinitesimal generator is well defined for 
every measurable and bounded function f : N* —> M by 

£f(n) = (f(n + 1) - f{n))bn + ( f(n - 1) - f{n))n{d + cn + Bh* n ) l n > 2 . (8) 

The proof of the convergence, stated in Section [4] relies on a compactness-identification 
technique. The result proven in Section [4j Theorem |4.2[ also states the convergence in law of 
the sequence of occupation measures associated with the fast component H £ (see [31]). The 
main interest of this result is that in the limit, the study of the prey-predator community is 
simplified, since it is entirely described by the one dimensional birth and death process N. 
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2.3 Pathwise construction and first properties 

Following Fournier and Meleard ([28]) we construct a trajectory of the prey-predator process 
Z as a solution of a stochastic differential equation driven by Poisson point processes. Let 
Qi(ds,du) and Q 2 (ds,du) be two independent Poisson point measures on M + x M + with 
intensity dsdu the product of Lebesgue measures. We define for any initial condition (n, h) £ E 
the coupled dynamics 


N t = n + 


l«<ft 7 V s _ Q i {ds, du) 



lu<N 3 -(d+cN s -+BH a -)l N > 2 Q 2 (ds, du), 


(9) 


j f Ht = H t ( rBN t -D- CHt ), H 0 = h. 


A unique solution of these equations exists as long as the number of individuals remains finite. 


Theorem 2.4. Under the assumption that there exists p > 1 such that 

E^(iVo) p ^ < oo and < oo, 


i) For all T > 0 

e(" sup ( N t ) p 'j < oo and E^ sup < oo, 

hs[o,t] ' Ve[o,T] ' 

ii) If p > 1, there exists a unique solution ( Z t )t >o £ D(M + ,£^) of |9]). Its infinitesimal 
generator is given by ([6]) for any bounded measurable functions f with /(n,.) £ C^(M + ) 
for all n £ N*. 

Moreover, the process Zt is a Feller process in the sense that for any g : E —>• M continuous 
and bounded, the function z i —> E z (g(Zt)) is continuous and bounded on E, Vt > 0. 

Hi) If p > 2, for all bounded measurable functions f with f(n ,.) £ C^(M + ) for all n £ N* and 
for all z £ N* x M + , 

M{ = f(Z t ) - f(z) - [ Af(Z s )ds, 

Jo 

is a L 2 —martingale starting at 0 with quadratic variation 



f(N s + 1, H s ) — f(N s , iL s )) 


2 

bN s 


+ 


(f(N a -l,H s )-f(N a ,H a )) 


N s (d + cN s + BH s )ds. 


Proof, (i) Let us remark that the process (Nt,t > 0) is stochastically dominated by a pure 
birth process (A 7 , t > 0) that jumps from n to n +1 at rate bn. From Theorem 3.1 in [28] we 
know that for all T > 0 E(sup fg [ 0 T ](W) p ) < 00 and thus E(sup t £[o,T](Nt) p ) < 00 . 
Concerning the predator density, we notice that for all (n, h) £ N* x M + the solution cj> n (h,t ) 
of ([3]) satisfies 

Vt > 0, 0 < (j) n (h,t) < h V /i*, (10) 
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Since ft* < rBn/C , we obtain that for all t > 0 


rB 

H t < Hq\/ — sup N s . 


C 


se[o,t] 


Then 

E( sup (. H t f) < E((ff 0 ) p ) + sup (N t y) < oo. 

\e[o,T] j ^ He[o,T] ' 

The fact that the infinitesimal generator is given by © and the proof of (in) can be easily 
adapted from [28] . 

It remains to prove that Zt is a Feller process. We adapt the method introduced by Davis 
[23]. The prey-predator community process differs from Davis’ setting since the jump rates of 
the prey population are not bounded. However, we overcome this difficulty using the moment 
properties given in (i). We denote by (Tj,T 2 ,...) the sequence of jump times of the prey 
population. It is always well defined since the jump rate admits a positive lower bound ft. Let 
g E Cb(E) and i/j E Cb(E x R+). We define the application G ^ on E x M + by 

G^(z,t) = E z {g(Z t )l t < Tl + y{ z T!,t - Ti)l t> T^j ■ 

Let ei = (1,0) be the first vector of the canonical basis on E and let us define a function 0 
on Ex M + by Q(z,t ) = /q 9(cf>(z, s))ds with 


9(z) = 9(n, ft,) = n(b + d + cn + Bh). (11) 

The function t 1 —>• 1 — e -0 ! 2 ’*) is the cumulative distribution function of the first jump time 
Ti conditionally on {Zq = z}. Then 


G^(z,t) = e e( - z,t ^g((f)(z,t)) + / / e 0(2,t) tl>(<f>(z, s) + e 1} t - s)bn+ 


0 JE 


s) — ei, t — s)n(d + cn + Bcf> n (h , s)) 


ds. 


Let us remark that z i->- G^(z,t) is continuous since z ^ 4>(z,t ) is continuous by Cauchy 
Lipschitz theorem for all t > 0, and the integrand is locally bounded. 

We now iterate the kernel G^. From Lemma (27.3) in [23] we get that \/k E N, 

G^izy) = G%^(z,t) =E z (^g(Z t )l t <T k+1 +y(Z Tk+1 ,t — Tk+i)lt >Tk+1 ^j ■ 

Then 

\Gf(z,t) — E z {g{Zt))\ < E z {J^\g(Z t )\ + \i/j(Z Tk ,t — 7fe)|)l t> T fe ) 

< (Iblloo + ||^l|oo)P*(Tfc < t). 

We deduce from (i) with p = 1, that the sequence of jump times (Tk)km converges almost 
surely to 00 , hence P 2 (7fc < t) — >k^-oo 0- To obtain the continuity of z 1 —> Ptg(z) it is 
sufficient to prove that the probability P z (Tfc < t) converges to 0 uniformly on compact sets 
of E. 

Let K be a compact set of E , and set N + = sup{n; (n, ft) E K} and H + = supjft; (n, ft) E K}. 
We construct a sequence of jump times (S)fc)fceN that stochastically dominates the sequence 
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of jump times (Tfc)fcgp$ for any initial condition in K. We start by bounding from above the 
prey and the predator populations. Similarly as above, we define a prey pure jump process 
(X t ) t >o starting from N + and a deterministic predator population process Yt starting from 


H+: 


X t = N + + 



^-u<bX s —Qi (ds , du) 


JYt = Yt (rBXt - D-CYt ), Y 0 = H+. 


The difference with point (i) lies in the fact that this coupling bounds from above every 
trajectory with initial condition in K : more precisely for any initial condition z E K, Nt < Xt 
and Ht <Yf for all t > 0, almost surely. 

We introduce a Poisson point process with intensity 6(X t ,Y t )dt and denote by (S'j)jgN its 
sequence of jump times. Since the rate function 9 increases, we deduce that for all z E E and 
t > 0, 

P z{Tk < t) < P(S fc < t). 


The probability P (Sk < t) converges toward 0 as k —> oo, since for all T > 0, E(sup sg [ 0jT ] X s ) < 
oo and E(sup sg [ 0j T] ^s) < oo. □ 


2.4 Derivation from an individual-based model 

In this part, we justify that the model 0-© derives from a microscopic model for the prey- 
predator community. We introduce a scaling parameter I\ tending to oo and consider that 
the number of predator is of order K while the prey number remains of order 1. At each time 
t >0, the microscopic community is represented by a vector (X f K , H where N f h ' E N is the 
prey number and e N is the number of predators. This process is a two-types continuous 
time Markov process whose transition rates are given for all (n, h) E N 2 by 


(n, h) —> 


(n + l,h) at rate nb 

(n — l,h) at rate n(d + cn + B K h) l n >2 

(n,h + 1) at rate hr K B K n 

(n, h — 1) at rate h(D K + C K h ). 


The parameters B K . r K , D K and C K are chosen as follows: 


B k 


B 

K' r 


K 


Kr, D k = D, C K 


C_ 

K' 


The predation and the competition among predators are normalized following [281 H2]. The 
parameter of conversion efficiency r K is scaled in order to maintain constant the benefit from 
predation. 

We consider the limit as K —> oo of the rescaled process (N K , ). 

Theorem 2.5. Assuming that the sequence of initial conditions 
{No, -ff- )k>o satisfies 

T E (K+f) 3 ) 
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< oo. 


( 12 ) 



and converges in law toward (no, ho) E E, then for all T > 0 the process (N K , ^~)k>o 
converges in law in D([0, T], E) toward the piecewise deterministic process (N t , H t ) defined by 
0-0 with initial condition (no, ho). 


The proof of this theorem is based on a compactness-uniqueness argument which derives 
from Theorem 3.1 in [19] and will not be developed here. The moment assumptions (12) 
ensure that the processes Z K and Z are well defined and that the assumptions of Theorem 
3.1 in [19] are satisfied. In the latter, the authors prove a similar result for a gene regulatory 
network in which the chemical reactions occur at slow or fast speed. 


3 Ergodic properties 

In this section, we study the ergodic properties of the prey-predator community process Z. 
We will prove the irreducibility of the process and of specific sampled chains. From these 
properties and a Foster-Lyapunov criterion, we will show that there exists a unique invariant 
probability measure and that the process is exponentially ergodic. 

3.1 Some definitions and known results 

Let us first recall some definitions. Let (Xf)t> o be a Feller process taking values in E a locally 
compact and separable metric space. We denote by C its infinitesimal generator and by Pt 
its semi-group. For every A E 13(E) we set ta = inf{f > 0,X t E A}. 

The process Xt is irreducible if there exists a a—finite measure v on E, called irreducibility 
measure, such that for all A E 13(E) 

u(A) > 0 ==> Vx E E, E,(/ lA(Xt)dtj > 0. 

The process Xt is Harris recurrent if there exists a a— finite measure \i on E such that 
VA E 13(E) 

p(A) > 0 => Vx E E, F x (ta < oo) = 1. 

A Harris recurrent Markov process is always irreducible (see }38j). 

Moreover, a Harris recurrent process Xt has an invariant measure n (see ED- In the case 
where this measure is finite, we say that Xt is positive Harris recurrent. 

For continuous time processes, the positive Harris recurrence can be derived from a Foster- 
Lyapunov inequality satisfied by the infinitesimal generator on some petite set. Recall that 
a set C C E is petite if there exist a probability measure a on M+, and a non degenerate 
measure v a on E such that for any z E C 

Pt(z, -)a(dt) > i'a(-). 

For an irreducible Feller process whose irreducibility measure has a support with non empty 
interior, all compact sets of E are petite sets (from Theorem 5.1 and 7.1 in [44] 1. 

We recall sufficient conditions for the positive Harris recurrence of a Feller process. 

Theorem A. (Theorem f.2 in / 39‘j) Let (Xf)t> o be a Feller process taking values in E. If 
the following conditions are satisfied: 
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i) X is irreducible with respect to some measure whose support has non empty interior. 

ii) Foster-Lyapunov inequality: there exist a function V : E —> [1, oo[ such that lim^i.^ V ( z ) 
oo, a compact set K C E and two constants 6 ,7 > 0 such that 

CV(z)<-^V{z) + 51 k {z), Vz € E. 

Then X is positive Harris recurrent and there exists a unique invariant probability measure 
n. Moreover tt(V) < 00 . 

The process Xt is ergodic if it has a unique invariant probability measure it and if 

lim | \Pt(x ,.) — tt\\tv = 0, Vx E E. 

t—> OO 

Moreover, Xt is exponentially ergodic if there exist a function R : E —> (0, 00 ) and 0 < p < 1 
such that 

|| Pt(x ,.) — 7r|| tv < R[x)p l , Vx € E. 

In the case of continuous time Markov processes on continuous state spaces, the ergodicity 
is related to the behavior of skeletons of the process. A skeleton corresponds to a sampling 
of the continuous time process at some fixed time. For all A > 0, the A—skeleton of X is 
the Markov chain (-XfcA)fceN with transition kernel Pa- We recall sufficient conditions for 
exponential ergodicity of a Feller process. 

Theorem B. (Theorem 6.1 in and Theorem 6.1 in \‘J9f ) Let (Xt)t >0 be a Feller process 
taking values in E which satisfies both conditions i) and ii) in Theorem [Aj If furthermore 
there exists an irreducible skeleton (A^A)fceN (^A > 0), then X is exponentially ergodic and 
there exist 0 < p < 1 and 0 < R < 00 such that, for all z E E, 

\\P t (z,.)-ir\\ TV <Rp t V(z). 

Let us briefly explain the origin of the condition on the skeleton of the process. Since Xt 
is positive Harris recurrent, the irreducible skeleton (X^a )fceN h as an invariant probability 
measure. Hence, the skeleton chain is positive recurrent and aperiodic (see Theorem 5.1 in 
P81). The irreducibility is crucial to obtain the aperiodicity. 

Moreover, from the Foster-Lyapunov inequality ii) in Theorem [Aj we deduce that the skeleton 
chain also satisfies a Foster-Lyapunov inequality with the same function V: there exist 7 ' < 1 
and 5' > 0 such that for every initial condition z G E 

E 2 (H(X A )) <^V{z) + S'. 

From Theorem 6.3 in m, we deduce that the skeleton (AfcA)fceN is geometrically ergodic. 
There exist 0 < p < 1 and 0 < R < 00 such that for all z G E 

\\PkA(z,.)-n\\ TV <Rp k V(z). 

The exponential ergodicity of the continuous time process then derives from the semi-group 
property. 
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3.2 Irreducibly 

In this section we study the irreducibility of (Zt)t>o in E = N* X M_|_. Let us highlight that a 
Borel set A G B(E ) can always be written as 

A= \J{k}xA k , 

k> 1 

where A k G £>(R_|_). We introduce the measure a on E as the product of the counting measure 
on N* and the Lebesgue measure A on R + : 

VAet3(E ), a(A) = Y^KA k ). (13) 

k> 1 

In particular, if cr(A) > 0, then there exist k G N* such that A (A k ) > 0. 


Theorem 3.1. (i) IfrB — D < 0, then the process ( Zt,t > 0) is irreducible for the measure 

a on E given by (13). 

(ii) Otherwise, the process (Z t , t > 0) is irreducible for the measure a' which is the restriction 
of a to the space E' = N* x [h\,+ oo), for h\ defined in ©■ 

The dichotomy in this result derives from the fact that when h\ > 0, the set N* x (0, h *) 
is transient for the dynamics. Indeed, the flows t H > fi n (h,t) are increasing functions for 
(n, h) G N* x (0, hfi) and therefore the trajectories cannot enter this area. 

In the sequel we prove a stronger result on the probability for the process Z% to reach open 
Borel sets, from which Theorem 3.1 follows. 

Theorem 3.2. (i) In the case where rB—D < 0, we consider an interval I = {k}x(h-, h + ) 

with 0 < /i_ < h- 1_ and k G N*. Then for every initial condition (n, h ) G E, there exists 
to > 0 such that Vi > to, 

{%t £ I) > 0 . 

(ii) We have a similar result in the case where rB — D > 0 for any interval I C E' such 
that a'(I) > 0 and any initial condition ( n,h ) G E'. 


The proof derives from the construction of ideal trajectories and from comparisons between 
the different predator flows. 


Proof of Theorem 3.2. 


(i) We assume that rB 


D < 0 which is equivalent to h\ = 0 . 


We consider different cases depending on the position of the interval I = {fc} x [h-, h+) with 
respect to the line n i —> h* of the predator equilibria and on the initial condition (n, h ) G E. 
These cases are illustrated on Figures [l] to [7j On these Figures, the state space E is represented 
as the positive quadrant of M 2 separated by the line n >-)• h* n . The process ( Z t , t > 0) can only 
cross this line by a jump of the prey number. When the process is above this line, the predator 
density decreases, while it increases when the process is under this line. The intervals which 
cross this line (i.e. such that ht G (h_, h + )) will play a specific role in the proof since they 
are stable by the predator flow cj) k - 

We introduce additional notations: for all m G N*, and i,j/ G M+ we set r m (x,y) the time 
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needed for the flow <j> m to go from x to y. This time is well defined for x > y > h* m or 
x < y < h* m . In these cases, it satisfies 

4>m(x,r m (x,y )) = y. 

First case: the interval I is stable for the flow 4>k (be. h k £ (h-,h + ) if h* k > 0 or 
h_ = h* k = 0 otherwise) 

Our aim is to prove that for any t > to, F^ n ^(Z t G I) > 0 for some to > 0. The idea is to 
construct simple trajectories which enter the interval I and arise with positive probability. 
We split the reasoning into different sub-cases depending on the initial condition. We focus 
on initial conditions such that h > /i*. The other cases can be treated similarly by symmetry. 

A) If n < k and h- V /t* < h. 

We first consider the specific sub-case where n < k and h- < h* n < h < h + (see Figure [I]). 
In this setting, we are interested in trajectories with exactly k — n prey births. These trajec¬ 
tories reach the line {k} X K+. Furthermore, the number of predators remains in the interval 
[h* n ,h V h k \ C (h-,h + ). This property derives from the fact that the predator density de¬ 
creases as long as Ht > h* Nt and remains therefore smaller than h but greater than h* N > /i* 
since Nt > n. If the process jumps below the line then the predator density increases 

and remains bounded by h k . Thus, after k — n births events, the process reaches the interval 
I. 



Figure 1: (Case IT.) Different ideal trajectories for the specific sub-case where the initial condition 
in, h) satisfies n < k and /i_ < h* n < h < h+. The red line is the map m h^. 

Let us now prove that such trajectories occur with positive probability. We first compute the 
probability that the first jump is a birth : 

P (n,h) (first jump is a birth) = E( nift )(E(l first jump is a birth |?i)) 

/ bn \ 

(n,/l) \0(?z, 4> n (h,Ti))J 

bn 

~ 6{n, 4> n {h, h V h* n )y 
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where the total jump rate 9(n, h) defined in © increases in n and h. 

Then by induction, the probability that the k — n first jumps are births, is greater than 

bn ^ b(n + l ) ^ ^ b(k — 1) 

6(n,h\J h* n ) 9{n + 1, h V ti*) 9(k — 1, h V It*) 

Recall that the sequence of jump times of the prey population is denoted by (T m )meN- Let 
us fix t > 0. Using the lower bound /i* of the predator population size, we bound from below 
the probability that the k — n births happen before t by 

P< t\k ~ n births) > P (Voiss(t9(n, h^)) > k — n) > 0. 

where Voiss{t6(n , )) is a random variable with Poisson distribution of parameter t9(n , /i*). 

Finally, we request that no other jump occurs before t, then 

P(n, h ) (Tfc_ n+ i > t\Tk- n < t and k — n births) > exp(— t9(k, h+)) > 0 

Then the event {T*, > t and Tk-n < t and k — n births} has positive probability, and on 

this event Z% G I. 


Let us come back to the general case where n < k and h- V /i* < h. 

We will consider the trajectories which remain on {n} x M + until Ht reaches h + . Then, we 
will request that k — n births occur before the predator population size reaches . 
Therefore, we define r\ = r n (h, h + ) when it exists and set r\ = 0 otherwise. The first step 
is to require that Tj > r\. Since in this case h + > h* k > and h > h * n , the flow < j> n (h ,.) 
decreases and thus 

W(n,h)(Ti > n) = exp^- J 9(n, (j> n {h, s))ds^ > exp (-n9(n, h)) > 0. 

Then, we define r 2 = r n (h + ,h-) when it exists (i.e. if G (h-,h + )) and set r 2 = +oo 
otherwise. It is important to remark that r± + r 2 = r n (h,h-). The specific case considered 
above corresponds to rq = 0 and r 2 = +oo. 

We request that Tk- n < n + ^2 and that these k — n jumps are births. An easy adaptation 
of the previous result shows that this event has positive probability. 

Moreover, at time T^_ n , Hx k _ n G The upper bound Hx k _ n < h + derives from the 

same reasoning as above. The lower bound of the predator density comes from comparisons 
of the different flows. We denote by F/ the vector field associated with (pf. F[(y) = y(rBl — 
D - Cy), for l G NT 

The first birth occurs at time r\ < Tj < rq + r 2 and Ht 1 = (j> n {h, Ti) > /i_. The second jump 
happens at T 2 < rq + r 2 and H T2 = (pn+iiH ^, T 2 — T{) > h-. Since F n+] (y) > F n (y) for 
all y G M+ then <p n+ i(H Tl ,T 2 — Ti) > (p n (H Tl ,T 2 — T{) and thus Ht 2 > cp n (h,T 2 ). Then by 
iteration, we deduce that Hx k _ n > (j> n (h,Tk- n ) > <l>n{h , r i + r 2 ) > h-. 

We define the time to = n + r 2 lr- 2 <+oo- Let us now consider t > to and finally request that 
Tk-n- i-i > t. As in the previous case, we deduce that these trajectories occur with positive 
probability and satisfy Z t G I. 

B) If n < k and h* < h < h-. (See Figure [3]) 
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Figure 2: (Case 1 A.) An ideal trajectory. The vertical blue arrows represent the time needed by the 
flow to get from the tail to the head of the arrow. 


The challenge is to increase the predator density up to /r_. Let us fix a time s > 0. We consider 
trajectories which have exactly k — n jumps before s, which are births. Then using a similar 
reasoning to case A), we deduce that H s > <f> n (h, s). We define the time r\ = r k ((f> n (h, s), /i_). 
Therefore, for every t > to = s + ri, if no jump occurs on the time interval [s, t], then 

H t = (j>k(H s ,t - s) > (j>k(H s ,n), 

since t — s > r \. Moreover, (f>k(H s ,ri) > 4>k((f>n(h, s), r\) = h- as H s > (j) n (h,s). Thus, Z t G I. 



Figure 3: (Cases IB. and 1C.) Examples of ideal trajectories. 

C) If n > k and h > h*. (See Figure [3]) 

The reasoning is similar to the previous case, except that we aim at decreasing the preda¬ 
tor density. We consider trajectories which have exactly n — k deaths before s > 0. Then 
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H s < (f> n (h,s). We define the time r\ = r k {cf) n {h, s),h + ) when <fi n (h,s) > h + and set r\ = 0 
otherwise. For every t > to = s + r\, if no jump occurs on the time interval [s, t], then Zf E /. 


Second case: The interval I is below the line n i —> h* n { i.e. h* k > h + ) 

We will construct an auxiliary interval which is stable for the predator flow. Then, we will 
prove that starting from this interval, the process enters / in some finite time. 

We introduce the integer m = max{^ E N *,h* < h+}. Once again, we split the reasoning in 
three cases depending on the position of h* m with respect to (h_, h + ) and the positivity of h* m . 
A) If /i_ < < h + . 

We define the interval J = {m} x [ }l++ ^ hm ; /*,_) which is stable for the flow <j) m (see Figure Q. 
Then, from the first case, there exists to, such that Vi > to, P (Z t 6 J) > 0. 

We set r\ = r k ( h+ ~^ hm ,h+). Let us remark that the trajectories starting from {no, ho) € J 
such that exactly k — m births occur during n satisfy that Z ri E /. This derives once again 
from comparisons of the flows <j>i for m < l < k. Moreover, such trajectories arise with positive 
probability. 

Therefore for any t > to + r\, we deduce from the Markov property at time f — r\ that 

W(n,h){ Z t £ I) > E (n,h) (^Zt- ri ( Z ri G I)^Z t - ri eJ^(n,h){ Z t-ri ^ j) > 0. 



Figure 4: 


(Case 2 A.) Construction of the auxiliary interval and of an ideal trajectory. 


B) If 0 < h* m < /m < h+. 

For this configuration we use the invertibility of the flow </>£. We will construct an interval 
I' = {k} x (u,v) such that (f>k(u,so) = h- and (j>k{v,so) = h+ for some so > 0 and that 
furthermore satisfies that u < < v. 

To this aim, we fix e E (0 ,h* n ) and remark that for any k > 1 and 0 < £ < y < h < h* k , 
the equation (j>k{y,s ) = h is equivalent to y = ipk^is), where ipk,h{s) is the inverse image by 
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Figure 5: (Case 2 B.) Construction of the auxiliary interval V. 


the flow i fik(.,s) of the point h. We deduce from Q that the application ipk,h is defined from 
[0, r fc (e, h)] to [e, h] by 


gs(rBk—D) _ hC ^Qs(rBk-D) _ 


It is continuous and strictly decreasing on [0, r k (e, h)}. 

Furthermore, from the uniqueness of the flow we deduce that for any 
r fc (e, h) > s> 0 

ipk,h + {s) > ipk,h-(s). 

Therefore, there exists a time so > 0 such that the points v = i/’k,h + {so) and u = ( s o) 
satisfy u < h* m < v and we set I' = {k} x (u, v ). 

From the case 2.A) we deduce that there exists t\ such that Vt > t\, 

P {n , h) {z t ei') >0. 

For any trajectory which is in I' at time t, we request that no jump occurs during sq, which 
happens with positive probability. Therefore, using the Markov property at time t, we deduce 
that 

^(n,h){Zt+s o C I) > 0. 

C) If h* m = 0 

In this case the above construction 2B) does not work because the only stable interval on 
{m} x M_|_ are of the form {m} x (0, a) with a > 0, which would impose so = +oo. 

Let us fix a small 5 > 0 and define the interval / 3 = {m} x (h_ + 5,h+ — 5). We remark that 
we can adapt the previous reasoning to prove that for every {n,h) G E there exists to > 0, 
such that Vt > to i 

p {n, h )(Zt e 1 3 ) > 0. 
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Let us explain this construction (see Figure ph. We first define m! = min{q E W,h* > 
h- +6}. As in the step 2B), we construct an auxiliary interval I 2 = {m} x (a, b ) with so > 0 

/l_|_ 5 — 4>m(b, So), 

h— T & — 0m(®, So)) 
and a < h* m i < b. 

We fix r) > 0 such that h* n , > a + r/ and set / 1 = { m'} x (a + rj,b). From the first step 
there exists t\ such that Vf > t±, IP ( n ,h)(Zt £ I 1 ) > 0. Starting from / 1 we request fur¬ 
thermore that exactly m! — m successive deaths occur on the time interval [t, t + n] with 
n = r m (a + 7?, a). This ensures that Z t+ri E I 2 . Furthermore we request that no jump occurs 
on [t + ri,t + ri + s 0 ], and thus, Vt > h V( n ^(Z t+ri+S0 E I 3 ) > 0. 

We now define the times r 2 = r m (h-,h- + 5) and r 3 = r k {h + — d,h+) and set £2 = 
min(?’ 2 , 7 - 3 ). For any trajectory which is at time t + r\ + So in I 3 we request furthermore that 
exactly k — m successive births occur before the time t + r\ + so + £2 • Therefore, we deduce 
from the Markov property that Vt > to 

^(n,h)(Zt+r\+so+t 2 G L) > 0. 

Third case: < h-. 

The proof is very similar to the second case. We introduce the smallest integer m such that 
h* m > and adapt the previous reasoning by inverting birth and death events. 

(ii) Let us now consider the situation where h\ > 0. Starting from a point (no, ho) E E 
such that h > h\, the process cannot reach the set {z E E,h < h\} which corresponds to 
the hatched zone on Figure [7} Therefore, we restrict ourselves to the measure a' and initial 
conditions in E 1 . The proof is the similar to above. 

□ 
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Figure 7: The hatched area corresponds to the points that can not be reached starting from an initial 
condition which is not in this zone. 


Proof of Theorem 3.1. We give the proof in the case where h\ = 0, the other case being an 
easy adaptation. 

For any A £ 13(E) such that a(A) > 0, there exist an integer k E N* and a Borel set A k E B(M) 
such that {k} x A k C A and A (A k ) > 0. Once again, we split the proof in two sub-cases. 
First case: Let us first assume that there exists an open interval (/i_, h+) C A k with h- < h+ 
and define / = {k} x (h-,h + ). 

We choose e small enough such that the interval L = {k} x (/i_ + e, h + — e) still satisfies 
o(I e ) > 0 and fix a small 5 > 0 such that V(k, h!) E I £ , 4>k(h',5) E I. 

From Theorem 3.2 we can construct trajectories that belong to the interval I £ at time t > to 
with positive probability for some to > 0. We ask furthermore that no jump occurs during a 
time 6. Then for any t > to 


P(n,ft) (Z s £ 3, Vs £ [t, t + 5]) 

> P t n ,h) (Zt £ E and no jump occurs on [t, t + A]) 

> e - se ^ h +-^F {n>h) (Z t £ I e ) > 0. 


Therefore ^ 

^(n,h) ^ J 1 A(Z s )ds) S j > (5P( raj /,) ( Z s £ I,Vs £ [i, t + $]) >0. 

Second case: We now consider the case where Ak doesn’t contain any interval (as an 
example M \ Q, or a fat Cantor set). We consider an open bounded interval (h_,h+) such 
that A(7L n (h_, h+)) > 0. Such an interval always exists since 

OO 

X(A k ) = J2 A (A k n (M, M + 1)) > 0. 

M =0 

Moreover, it is possible to choose (h_,/i + ) such that h £ ^ [h-,h + \, i.e. this interval is not 
stable for the flow <f>k- Indeed, the opposite case would imply by successive divisions of 
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the interval, that for any e > 0, A (A k \ (h* k — e, h* k + e)) = 0. Thus, we would have that 
0 < X(Ak) = A (Afc n (h* k — e, h* k + e)) < 2e for any e > 0 , which is not possible. 

We now restrict ourselves to the set B k = A k n (/i_, h + ) with A (B k ) > 0 and assume that the 
flow 4>k increases on (/&_, h+) (the other case being an easy adaptation). 

Let us fix e > 0. In the sequel we consider the trajectories that reach the interval ( h _ — e, /t_) 
and then, we ask that no jump occurs until these trajectories attain h + . Then, the time spent 
by those trajectories in B k will be positive since the flows are continuous. More precisely, from 
Theorem 3.2 we deduce that there exists to such that Vt > to, T(n.,/i) {Zt £ {k} x ( h _ —e, /t_)) > 
0. We define the positive time r\ = r k (h- — e, h+) needed for the flow (ftk to go from /i_ — e 
to h- 1 _. For all t > to, we consider the event 

£t = ^Zt E {A:} x (/i_ — e,h-) and no jump occurs on [t, t + n]|. 

Then (ft) > 0 and 

(n,h)(^J l{k}xB k {Z s )ds \ St'j = ^(n,h)(^J 1 B k ( ( t ) k{Zti s))d. 


E 


:a I ft 


Since the flow cj) k is invertible, we make the change of variable u = (f> k (Zf . s ) and obtain 

rt+ri 


a t+ri 

^-{k}xB k (Zs)ds | f t 


= E 


(n,h ) 


MZt,n) 

z t u(rBk — D — Cu ) 




du | ft 


Since, Zt e (ti_ — e, /t_) and 4>k(Z t , ri) E (/i+, /i^), we deduce that for some constant v > 0, 

^(n,h)(J 1{ k}xB k (Zs)ds | ft) > vX(B k ) > 0. 

Thus, E (ni/l) (/“ ls fc (^s)ds)) > 0 which concludes the proof of the irreducibility. 


□ 


3.3 Positive Harris recurrence 

We recall the expression of the infinitesimal generator of the prey predator process given in 
©• We prove in the following that it satisfies a Foster-Lyapunov criterion. 

Proposition 3.3. Let V : E i —> [l,+oo[ be the function V(n,h ) = n + h/r.Then there exist 
5, 7 > 0 and a compact set K such that 

AV(z) < - 7 V(z) + 81 K (z), \/zeE (14) 


We combine Theorem 3.1 and Proposition |3.3[ to deduce from Theorem [A] in Section 3.1 that 


Theorem 3.4. The process Zt is positive Harris recurrent and thus there exists a unique 
invariant probability measure ir on E' which furthermore satisfies vr(P) < oo. 
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Proof of Proposition\3. S\ For all (n, h) G E 


AV (to, h) = —(rBn — D — Ch ) + bn — n(d + cn + Bh)l n >2 
r ~ 

C D 


= 


— h(h - \ — nine — (b — d)), if to > 2 

\ r r / 

— h(h -s') +6, if n = 1 

V r r ) 

For any 7 > 0, we obtain easily that there exists no G N* and ho > 0 such that 

— n(nc — (b — d)) < — 7 to, Vn > too 

and — h(h -) < — 7 —, V/t > Tiq- 

V r r J r 

Then for all (to, h) with to > no and h> ho, AV (n, h) < —yE(n, /<,). Moreover since no is finite, 
there exists h\(no) such that V(n, h) with n < no and h > /ri, AV(n,h ) < — 7 V(n,h). Simi¬ 
larly, there exists n\(ho) such that V(n, h) with h < ho and n > ni, AV(n, h) < —7 V(to, h). 
Therefore, setting A" = {1, • • • , max(no, ni)} x [0, max(/io, hi)} and 

6 = rna x{V(n,h),(n,h) G A"}, the function V introduced above clearly satisfies the Foster- 
Lyapunov criterion (14). 0 


Remark 3.5. The function V introduced in Proposition \3.3\ is not the only Lyapunov function 
of the system. An easy computation leads to the fact that W(n, h ) = n 2 + h also satisfies (14). 


With Theorem 3-4 we conclude that 7 r(VF) < 00 which improves our knowledge of the invariant 


probability measure it. 


3.4 Exponential ergodicity 

In this section we investigate the convergence in total variation norm of the transition kernel 
toward the invariant measure. Let us first recall Theorem 12.21 


Theorem. [Theorem 2.21 The community process (Z t )t >0 is exponentially ergodic. It con¬ 


verges toward its invariant probability measure 7 r at an exponential rate. There exist 0 < p < 1 
and 0 < R < 00 such that, for all z G E', 


| Pt{z, •) - tt\\tv < Rp t V(z ) 


(15) 


Proof. From Theorem |B] in Section 3.1, it remains to prove that a skeleton chain of the prey- 
predator process is irreducible. This condition is actually equivalent to the ergodicity for 
positive Harris recurrent processes (Theorem 6.1 j38j). 


It derives immediately from Theorem 3.2, that any skeleton chain (ZfcA)fceN (with A > 0) 
reaches any open Borel set with positive probability. Indeed, let O be an open set of B(E') 
with A (O) > 0 then there exists an interval / C O with A (I) > 0 and thus, for any z G E' 
and for q large enough, f’ z (Z q/ \ GO)> P z (Z q a G I) > 0. 

To generalize from open Borel sets to Borel sets, we need some regularity of the function z e->• 
1 F > z(-^a G A) for A G B(E’). We compute this probability by distinguishing the trajectories 
with the number of jumps J(A) occurring on [0, A], then for all z G E' and A G B(E'), 
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( 16 ) 


P z (Z A £ A) = J2 f z (Z a £ A and J(A) = k). 

k =0 

We recall that the sequence of jump times of the prey population is denoted by ( Tk)k£N and 
that P z (Ti > t) = where 0(z,t) = d((f>(z, s))ds and the total jump rate 9(z) is 

given by (11). 

The first term of ( |16[ ) handles trajectories where no jump occurs. It is given by 
F z (Z a g a and J(A) = 0) = e~ e ^ ds l A ((/>(z, A)). 


This function is not continuous in z since the indicator function 1 A is not continuous and the 
total flow cj)(z, A) is continuous. 

The idea is then to bound from below P z (Z A G A) by a continuous function (see Chapter 6 
of [KOj and 0)- In the sequel we consider 

¥ z (Z A G A) > P z (Z A £ A and J( A) = l), 

and prove that for any A £ B(E'), the function z i-)- T(z,A ) = P 2 (Za £ A and J(A) = l) is 
continuous on E'. The continuity will derive from the fact that the law of first jump time has 
a density with respect to the Lebesgue measure. 

Indeed, for any z = ( n , h) £ E’ 

^z(Z A £ A and J(A) = l) = P 2 (Z A £ A and Ti < A < T 2 ) 

A , N 

l A [n + l,c/)n + i{Mh,s), A - a)J bne -Q(.n,h,s) e e(n+l,Mh,s)A-s) ds 

+ J 1a (n - 1 , (t>n-i (4>n(h, s), A - s)^jn{d + cn + B4> n (h, s))l n > 2 
X e -9(w,ft,a) e &(n-l,<t>„(h,s)As) fig 



The first integral corresponds to the event where a birth occurs at T\ while the second interval 
to the event where a death happens at T\. In the sequel, we consider the first integral. The 
study on the second integral is very similar and will not be detailed. The predator density at 
time A conditioned on the fact that only one jump happens on [0, A] and is a birth occurring 
at time s £ [0, A] is given by 


9(n,h, A)(^0 — (ftn+l s), A s) . 


We note that for any s and A, the application (n, h) i-)- g^ n ,h, A)( s ) is continuously differen¬ 
tiable. To perform the change of variable y = g^h, A)( s ) in the previous integral, we have to 
verify that ^ 9(n,h,A )( s ) does not vanish. 


d_ 

ds 


9(n,h,A) (s) = d 2 4>n(h,s) ■ di4> n+ i((/) n (h, s), A 


- s) - d 2 (/>n+ 1 (' <t>n{h , s), A 



We recall that 4> n is the flow associated with ([ 3 ]), then 


d 2 (j> n {h, s) = (j) n {h, s)(rBn - D - Ccf> n (h, s)). 
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From the exact expression © we obtain that 


di cj>n(h,s) = 


(,1>n{h,S ) 


K l + J^{z {rBn - D)s - l )Y 

Then an easy calculation using Q leads to 


j~ s 9(n,h,A)( s ) =-rB 


9(n,h, A) ( , ®) 


1 i <Pn(h,s)C / (rB(n+l)-D)(A-s) _ ]\ 
' rB(n+l)—D \ / 


< 0. 


Let us finally remark that 


9(n,h, A)(°) = 4>n+i{h, A) and 5(n,fc,A)( A ) = <l>n(h, A), 


hence 


j\ A (n + l, cf n+1 (^ n (/i, a), A - a)) 6ne - e (-^^) e e («+h^(^ S ).A-,) ds 


r<j>n+i(h, A) 

/ lA(n + l,y)/(n,/i, A,y)dy, 

Jd)n.(h,A) ' ' 


(17) 


'4>n{h, A) 


where 


f(n, h, A, y) = frra- 


= -0 (n,M {n ( hiA) (3/)) e (n+l,0„ (M (n ^ A) (3/)) A~9 (r f thtA) («)) 


ds9(n,h, A)( s ) 


~ 9 (n,h,A)(y) 


The function of (n, h) given dehned by © is continuous since the upper and the lower bounds 
of the integral are continuous functions of (n, /i) on F7 and the integrand is continuous in (n, /i) 
on E' and locally bounded. 

To conclude with the irreducibility of (.ZfcA)fceN> we fix a point zq G if' and remark that the 
measure T(zo, •) = P z (Za £ ■ and J(A) = l) is non degenerate since T(zq. E') > 0. 

For any A G B(E') such that T(zq, A) > 0, there exists, by continuity of T, an open neighbor¬ 


hood O of zq such that Vz G O, T(z, A) > T(zq, A)/2. Moreover we deduce from Theorem 3.2 


that for any initial condition (n, h) there exists q G N such that P( n ,/i)(^gA G O) > 0. Then it 
derives from the Markov property at time qA and from the properties of the kernel T that 

^(n,h){ Z (q+ 1)A e -4) > ^(n,h){lz iq+1)A £A^Z qA {lz A £o)) 

> p m ( VgO )^>0 


Hence, (Zf.A)keN is irreducible with respect to the measure T(zq, •). 


□ 


Remark 3.6. Here we chose to study the behavior of skeletons of the process in order to 
derive the exponential ergodicity of the process. This standard method has already been used 
in the context of PDMP (see for example mw- Costa and Dufour developed a 

different approach based on an embedded Markov chain whose long time behavior is equivalent 
to those of the PDMP. For the prey-predator process Z it seems to us more natural to study 
the continuous time trajectories since comparisons where possible. Moreover in future work, 
we aim at deriving from these trajectory constructions more quantitative convergence speeds 
using recently developed coupling methods for PDMP (e.g. Mm- 
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Remark 3.7. The results of this section can he extended to the model including migration 
introduced in Remark \2.1\ Indeed, using the very same trajectory construction, one can prove 
that the process (j\d m ), H( m ' ) ) is irreducible on Efor the Lehesgue measure a. Moreover, 
the infinitesimal generator (J7J) satisfies a Foster-Lyapunov inequality with the same Lyapunov 
function. Therefore (N^ m \ H^) admits a unique probability invariant measure. Finally, by 
adapting the proof of Theorem 2.2 we get that it is exponentially ergodic. 


4 Re-scaling the predator dynamics 


We introduce a new parameter e which rescales the predator dynamics and illustrates the 
biological assumption that the predator mass is almost negligible comparing to the prey 
mass. There exists an important literature about the metabolic theory which describes the 
relationships between mass and metabolic characteristics of living individuals (see among 
others eh nm). Following this theory, the demographic parameters of individuals increase 
when their mass decreases. Here we simplify these relationships by assuming that the predator 
parameters vary as 1/e. 

For any e € (0,1], we consider the community process Z £ = (N E , H £ ) where for all t > 0, 


d H £ 

_ H f = -L(rBN t £ -D-CHf) 


(18) 


and the dynamics of N £ is given by the jump mechanism ([2]) associated with the predator 
population H £ . The process studied in the previous sections corresponds to e = 1 or to the 
parameters r £ = r/e, D £ = D/e and C £ = C/e. 

This scaling changes the time scale of the predator flow. If is the flow associated with (18) 
then 

4>n(h,t) = (^(h, 1 -), Vt > 0, V(n, h) 6 E. (19) 


4.1 Convergence toward an averaged process 

In the sequel we study the limit as e tends to 0 of the sequence Z £ in D([0, T\,E). The 
prey-predator process is a slow-fast system. As £ diminishes, the predator process converges 
faster to its equilibrium between the jumps of the prey population. The slow dynamics of the 
prey population is then averaged on the predator equilibria. 

We first give the expression of the infinitesimal generator of the process Z £ defined above: 
V/ E £ and (n, h) G E 

A £ f(n,h ) =d 2 f(n,h) h ^ rBn ® ^ + bn(f(n + 1 ,h) - f{n,h )) 

+ n(d + cn + Bh)(f(n - 1 ,h) - /(n, h))l n > 2 . 

To carry out the limit as e —> 0, we use the fact that the convergence speed of the flow 
<f> £ n to its equilibrium /i* is uniform in e. This is only true if the number of predators remains 
bounded from below by some strictly positive constant. To this aim, we make the following 
assumptions 
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( 20 ) 


i) rB — D > 0. This implies that the predator population cannot 
become extinct 

ii) We restrict ourselves to initial conditions in the set 

^ x [h{,oc). 

The state space E' is stable for the prey-predator dynamics and is the support of the irre- 
ducibility measure a' introduced in Theorem 3.1 ii). 


Proposition 4.1. Under Assumption (20), there exists 5 > 0 such that for all ( n,h ) E E', 

W n {h,t)-K\<\h-h* n \e- & -e 


( 21 ) 


Proof. Let us remark that thanks to (19), it is sufficient to prove the result for e = 1. 
Using an easy computation leads to 

h - ht 


|<MM) - K\ = 


1 + A( e -t(rBn-D) _ -Q 


, V(n, h ) 6 E' and t > 0. 


Therefore it is enough to find <5 > 0 satisfying V(n, h ) 6 E' and t > 0: 


e dt < 1 + E— ( g—t(rBn—D) _ 

K 


i). 


The function t i —> 1 + jfr(e t ( rBn D ) — 1) — e St equals 0 for t = 0, and increases as soon as 
hC > 5 and rBn — D > 5. Since h > h\ > 0 we choose <5 < Ch\ to obtain (21). □ 


Following Kurtz m, we introduce the predator occupation measure 

T £ ([0,t],^4) = [ 1 A {H e s )ds, Vt > 0 and A e B(R + ). 

Jo 


( 22 ) 


This random measure belongs to the set A4 m (M_|_) of measures p, on M + x M + such that 
/i([0,t] x M + ) = t, Vt > 0. For any t > 0, we denote by the set of the measures 

p E A4 m (M+) restricted to [0,f] X M + . 

In the sequel we prove using the averaging method developed in m that the sequence {N e , T £ ) 
converges in law. This method allows us to avoid the difficulties related to the fast convergence 
of the predator flow to its equilibrium. 


Theorem 4.2. FixT > 0 and assume (20). We suppose that the sequence of initial conditions 
(Z £ o) o<e<i converges to Z o in law and moreover that 


sup E((JVq) 4 ) < oo, sup E((FIq) 4 ) < oo. (23) 

0<e<l 0<e<l 


Then the sequence [N £ , T e ) converges in law toward ( N,T ) in B([0,T],N) x A4^)(M+) as 

e 0. _ 

The process N is a pure jump process on N* whose infinitesimal generator is well defined for 
every measurable and bounded function f : N* —> M by 


£f{n) = (f{n + 1) - f(n))bn + (/(n - 1) - f{n))n(d + cn + Bh* n )l n > 2 . (24) 

Moreover, the limiting measure is defined by T(ds X dy) = 5h*_ ( dy)ds. 

Ns 
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We say that N is an averaged process since it behaves as if the predator density is constant 
at its equilibrium. Let us consider the specific case where D = 0. In this case, the averaged 
prey population N evolves almost as a logistic birth and death process with individual birth 
rate b and individual death rate d + cn where c = c T rB 2 /C. The logistic parameter c 
corresponds to the apparent competition pressure (see my- it takes into account both the 
prey competition c and the effect of predation. 


Proof. We divide the proof in several steps. The first three steps are devoted to the conver¬ 
gence of the prey population process. We use a standard compactness-identification method 
inspired from Genadot [29]. In the fourth step, we use the averaging method developed by 
Kurtz m to prove the convergence of the predator occupation measures. 

Step 1: We prove that there exists an unique (in law) solution to the martingale problem 
associated to ( [24] ) : for every measurable and bounded function / : N* —> M 

f(N t ) - /(JVo) - f Cf{N s )ds 
Jo 


is a martingale. 

It derives from the representation Theorem 3.2 in [32] and a localization argument (see Theo¬ 
rem 4.6.3 in [26] ) that the uniqueness of the solution of this martingale problem is equivalent 
to the uniqueness in law of the solution N of the following stochastic differential equation 
provided that sup s<t N s < oo a.s. for any t > 0: 



N t - N 0 + , , ■ L u<bNa - 

1 0 JR+ 


L bN 


<u<N a -(b+d+cN a - +Bh— )Q{ds,du), 


(25) 


where Q is a Poisson point measure on (M + ) 2 with intensity the product of Lebesgue measure 
dsdu. The uniqueness of the weak solution of (25) can be adapted from [28l 

and sup s<t N, 


E(IVo) < Too, then this solution is well defined on 


Moreover, if 
< oo a.s. for any t > 0. 


Step 2: Tightness of (Nf, 0 < t < T) 

Similarly to the proof of Theorem |2.4[ we construct a pure birth process X independent from 


Hq and thus from e, which dominates the prey population. Then we deduce from (23) that 
Vz € E' 

sup E z ( sup (IVf) 4 ) < E-( sup (W,) 4 ) < Too. (26) 

se[o,T] se[o,T] 


0<e<l 


Let us now fix 77 ,<5 > 0 and consider stopping times a, r such that a < r < (a + 6) At. Using 
the trajectory’s construction 0 , we write 

Nf = Xf T [ Ng(b - (d T cNf + BH £ s )l N e> 2 )ds + Mf , 

Jo 

where Mf is a pure jump martingale with quadratic variation 


(M e ) t = / Nf(b+(d + cNf + BHf)l NI > 2 )ds. 
Jo 
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Hence 


e((J+-A+ 2 ) 

= E (( [ N £ (b - (d + cNI + BH £ s )l NS > 2 )ds + M £ T - M £ j‘ 


< 2 


E 


N £ s (b + d + cN £ s + BH £ s )dsf) + E( (M £ - M £ ) 


We deduce from (10) that the first term is bounded above by 

rB 


E 


V LJ q \ 

N £ (b + d + cNg + B(Hq + — sup N £ )ds) ). 

C o,^rn cl / 


wG[0,s] 


Therefore, using (23) and (|26|), there exists a constant Ct > 0 such that 

E( 


:((J N £ (b + d + cN £ + BH £ )ds ) 2 ) < C T S 2 


For the second term 


E {[M £ - M;) 2 ) = e(^ iV|(6 + (d + cN £ s + BH £ )l]ye> 2 )ds 

< E ( [ N £ (b + d + cN £ + B(H £ 0 + ^ sup N £ ))ds ), 

Vcr uGfO.sl ' 


where the last inequality derives from (10). Combining this two inequalities (23) and (26), 
we deduce that 


E ((N £ -N £ ) 2 ^) <5C t , 


for some constant Ct > 0 independent of e which leads to the tightness of the laws of 
(Nf,t G [0,T]) in B([0, T],N) using Aldous tightness criterion p]. 

Step 3: Identification of the limit 

Let us consider a sub-sequence of (!V t e ,0 < t < T) (still denoted N £ ) which converges in law 
when e —> 0 toward a process (N t , 0 <t<T). In the sequel we prove that N is the unique 
solution of the martingale problem associated to (24). 

We consider bounded functions f, g±,..., gk onW and times 0<t\<---<tk<t< t+s < T. 
We deduce from Theorem 12.41 that 


E 


+E 


rt+s 


N t £ +S ) - f{N £ ) - / Cf{N £ u )du) H 9i{N £ ti ) 

Jt i =1 

(f(N £ - 1) - f(N £ ))N £ B(H £ u - h* N e)l N e> 2 du) n^(^) 


t+s 


(27) 


= 0. 


2=1 


From the convergence in law of N e to N and (26), the first term converges as e —> 0 to 

k 


E 


rt+s 


(. f(N t+s ) - f(N t ) - / Cf(N u )du ) n 9i (N ti ) 

Jt i=i 
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Let us prove that 


lim E 

£—^0 


t-\-s 


(f(K - 1) - f(K))KB(Hl - h* N e)l N e> 2 du) n 9i {Nl) 


2—1 


= 0. 


(28) 


We first remark that 

rt-\-s 


E 


< 


{f(K - 1) - f(K))KB(Ht - h* N ') 1 Ng > 2 du 

%— 1 

k nt-\-s 

H H^llooE / | f(K - 1) - f{N £ u )\N £ u B\Hi - - |1 N *> 2 du 


2—1 




We split the integral according to the sequence of jump times (T|, k E N) of the prey popu¬ 
lation: 


E 


a t+s 

I f(K - 1) - f(K)\KB\HZ ~ h* N e\l N e> 2 du 
00 00 / r T k+ i A (*+ s ) 

<2||/||ooS EE 


,T £ 




fc=0 n=l 

Since on the event {N^e = n}, for any u E [T|, T|_ ) _ 1 A (i + s)], the predator density is given 
by H= (j) e n {Hj,e, u — T|), we derive from Proposition 


4.1 


that 


S(u-T^) 

\h* m -H £ u \<\h* Nf -HUe -— 


Hence 


E 


a t+s 

(WZ - 1 ) - f(K)W u B{Hi - h* N e)l N e> 2 d,U 
00 00 / r T Ui a(*+s) 

<2||/||oo5 EE E( l n (A^ 


fc=0 n=l 






S(u- T %) 

\h* n - Hfr\e - ~du 


From (10), we obtain that for every T > 0 


rB 


sup H u < H + + — sup 1\P < +oo. 
mS[ 0,T] ^ ue[0,T] 

Since /i* < rBn/C, we deduce that for all u E [T^, T| +1 A (t + s)] 

\h*N£ ~ H-tA — H+ + 2— sup . 
T fc ° ue[o,T] 

With a change of variable u —>■ u — TE , we obtain that 


(29) 


f 

■>n 


T k+ iA(t+s) s(u-ri) r T k+ i A 0+ s )- T fe 

e-/ exp(--)du 

4o e 

/■ t+s <5u 

< / exp(-)dul T f<t+ s 

4o e 

e 
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Combining these two inequalities, we have that 


E 


(f(K - 1 ) - fiK))N £ u B(Hl - h* N e)l N e> 2 du 

/ rB °° \ 

< 211/Hoo^ EUH+ + 2— sup N £ ) sup N £ ^ lr'<t+ a ) 
v 0 ue[o,T] «e[o,X] r—r 7 


It remains to bound the expectation independently of e. To this aim, we denote by the 
number of jumps before time t. By neglecting the non-positive terms, we deduce from the 
trajectorial construction ([ 9 ]) that 


sup NfJj* < f [ (N £ s _ + .J £ _ + l) l u < bN i_Qi(ds, du) 
te[o,T] Jo Jr+ 

+ / / lu<NI_(d+cNj;_+BHI_)Q2{ds, du ). 

JO J K+ 


Then, we choose S > T and taking expectations we deduce from (23) and (26) that there 
exists a positive constant Cs such that 


sup <Cs + Ef f bJf sup N £ dt 

ue[o,T] 7 V 0 se[o,t] 


We conclude using Gronwall lemma that E^sup tg [ 0 T ] N^Jj,J < 00 . A very similar computa¬ 
tion leads to 


E( sup ( N £ fj f) < E( [ T 2b(Nl) 2 Jl + bNgJg + WVf(1 + N £ s ) 
ue[o,T] 7 v 0 

+ ( N £ s - 1 ) 2 N £ s {d + cN £ s + BH £ s )ds ) . 


We conclude with (23), (26) and Gronwall lemma that 

E^ sup < 00 . 

v te[o,T] 7 

Therefore, we deduce that there exists a constant Cs.t > 0 such that 

rt-\-s 


E 


(J \f(K - 1) - f(N £ ))N £ B(H £ u - h* N e)l N e> 2 du) | < Cs,T 


Thus, (p8|) is verified. 

E 


Therefore 


rt-\-s 


N t+S ) - f(N t ) - / CfiNjdujUgiiNt.) 


i= 1 


= 0 , 


and thus A 7 is a solution of the martingale problem associated with C. From (23) we deduce 
that E(iVo) < 00 and then from the first step, N £ converge in law to the unique solution of 
the martingale problem associated with £ which is a birth and death process which jumps 
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from n —> n + 1 at rate bn and from n —> n — 1 at rate n(d + cn + Hh*) l n > 2 - 


Step 4: Limit behavior of the predator population size 

We prove the weak convergence of the sequence of occupation measures T £ introduced in (22). 
From, Lemma 1.3 in m , the sequence of the laws of T £ is tight in the set ,M^(K+) if for any 
6 > 0 and t E [0, T] there exists a compact K C M+ such that 

inf E(r £ ([0, t] x K)) > (1 -5)t. 


From (29), we deduce that the family {H £ ,e E (0,1], s E [0,T]} is relatively compact, and 
thus the sequence of the laws of (r £ ) eg ( 0j i] is tight. 

Hence, the pair (!V £ ,r £ ) is tight and we consider a sub-sequence, still denoted (JV £ ,r £ ), that 
converges toward (N,T). 

Since r £ ([0,f],M + ) = t. for all e E (0,1) and t E [0, T], we deduce (see Lemma 1.4 in 
EH) that there exists a process taking values in the set of probability measures on M + , 
measurable with respect to (w, s), such that for every measurable and bounded function h on 
[0, T] x M + 

/ / h(s,y)T(ds X dy) = / / h(s,y)'y e (dy)ds. 


Using (26), we deduce from Theorem 2.1 in m that for every function / on N* continuous 
and bounded and t < T 


f(N t )-f(N 0 ) - 


(N s + 1) - f(N s ))bN s 


(30) 


/0 JR+ 

+ ( f(N s - 1) - f(N 8 ))N 8 (d + cN s + By) 1^> 2 ] ls (dy)ds. 

is a martingale. 

Concerning the predator dynamics, for every function g E C^(i?+), the process 

eM £ g (t) = eg(H £ ) - eg(H £ 0 ) - f g (H £ )H £ (rBN £ -D- CH £ s )ds- 

Jo 

is a martingale. Using (26), we prove easily that (eM £ (t)) Q<£<1 tg [ Q is uniformly integrable 
and that 

limE(|eM|(t) — M g (t)\) = 0, 


£—^0 


where M g (t) = J R+ g'(y)y(rBN s — D—Cy) , y s (dy)ds. Therefore we deduce from the uniform 

integrability of sM £ that M g (t ) is a martingale. Since it is also a continuous and finite variation 
process, it must be null. Thus 


/ E+ 


g'(y)y(rBN t — D — Cy)'yt(dy) = 0, for dt — almost every t E [0,T]. (31) 


(30) we deduce that for all t E [0,T], 


We recall that the infinitesimal generator of N is given by (24), then by identification in 

Vlt(dy) = h* N 


(32) 
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Then we apply (31) to the function g(y) = J°° f(u)du , where / is continuous on a compact 


of) 




/ f{y)y{rBN t — D — Cy)jt(dy) = 0, for dt — almost every t £ [0, T], 

From Riesz Theorem, the measure jit(-) = /' y(rBNt — D — Cy)^t{dy) on M + is null for almost 
every t £ [0, T\. Then, for almost every t £ [0, T], 7 1 only charges 0 and h*^ . Finally using 
(32) we conclude that 7 t.(dy) = 5h*- ( dy ) for dt— almost every t £ [0, T]. 

— JV t 

□ 


Remark 4.3. T/ie proof of Theorem \4- ^ relies on the uniform convergence of the flows (j> n (h, t ) 
to their equilibrium h* n as t —»• oo, V(n,/i) £ E'. It is not possible to ext end this reasoning to 
the process (N^ m \ H^) including migration introduced in Remark 2.1, since its state space 
is N x M+. However, we can still prove using the same method that the accelerate d seq uence 
(jVM^ r (m) 

,£ ) is tight (with obvious notations). The difficulty to extend TheoremUHa to the 
migration case lies in the identification of the limiting values. 


4.2 Long time behavior of the averaged process 

We are interested in the long time behavior of the averaged prey population N. 


Proposition 4.4. The process N is positive recurrent on N* and converges toward its unique 
invariant probability measure jl = Tnbn which satisfies the system 


Vn > 2, 


b n ~i 

^ n n II” =2 (d + ci + Bh*)^ 1 


(33) 


The proof is very classical and left to the reader (see for example |2], p.216). Moreover, we 
are interested in the process Z = (JV, hE) which represents the averaged prey and predator 
populations. We obtain immediately the form of its invariant distribution since h* N is a 
function of N. 


Corollary 4.5. The process Z = (N, h* N ) admits a unique invariant probability measure ir 
given by 

71 = ^ ] Trid(n,h* ) i ( 34 ) 


n> 1 


where the sequence y n satisfies (33) 


The expression of W is not explicit, however thanks to (|33j) we can derive information on 
the shape of the distribution. In the sequel we denote by \x~\ = min{m £ Z;m > x} for 
x £ M. 


Proposition 4.6. The invariant probability measure If admits a unique maximum at 1 if 
y\ > y ,2 and at n\ > 1 otherwise. The value of n\ can be explicitly computed in function of 
the model parameters. 
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Proof. The study of the existence of maxima of the distribution n and thus fx is equivalent to 
we study of ^" +1 — 1 = £(n) for £ the function defined by 

^ ~~ (x + 1) (d + c(® + 1) + 5(r5(x + 1 )-D)/C) ~ VX G + °°' ) ' 

The sign of £(x) is given by the sign of the polynomial 

ax 2 + /5x + 7 = 0 (35) 


with 


Its discriminant equals 


a = c + 


rB 2 

”cT 


P 


d +2c+ 2 


rB 2 


BD 

~C~ 


'y = -b + d + c + 


rB 2 
C 


BD 


P 2 ~ 4«7 = (d - ^t?) 2 + 4(c + ^r)b, 


C 


C 


which is always positive. Therefore the polynomial (35) admits 2 real roots: 

-P- v 7 ( P 2 - 4cr y) , _ -P + vW-4ot) 


x 0 = 


2a 


and xi = 


2a 


The smallest root xo is always negative. If x\ > 1 then the invariant distribution JI admits 
exactly one mode at n\ = [xi]. Otherwise, the sequence (jx n )n> l is decreasing. To conclude 
the proof it remains to remark that the condition x\ > 1 is equivalent to ^(1) > 0 which is 
the condition given in the proposition. □ 


4.2.1 Numerics 

For each e E (0,1] we proved in Section [3] that there exists a unique invariant probability 
measure n £ for the process Z £ . In this section, we study with numerical simulations the 
behavior of the sequence of invariant probability measures (tt £ ) £ as e —> 0. An approximation 
of the invariant measure n £ is obtained by simulating 3000 times the prey-predator process 
Z £ on a long time interval. One interest of the process Z £ , is that the flow for the fast 
predator population admits an explicit formula Q. Therefore, we are able to compute exact 
simulations of the process Z £ for every e E (0,1]. The code for these simulations is available 
in [IT ] (Chapter 3, Appendix A). 

In the simulations, the demographic parameters are given by 


b = 0.4, d = 0, c = 0.005, B = 0.02, r = 2, D = 0 and C = 0.04. 


In Figure[8[ we draw three-dimensional histograms of the distributions of the invariant measure 
ir £ for different values of e (in (a) e = 1, in |(b)| g = 0.1 and in 
that the support of these measures concentrates as e 


(c) e = 0.00001). We observe 
0 on the set {(n, h*),n E N*}. This 
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(a) e = 1 


(b) s = 0.01 



(c) £ = 0.00001 


Figure 8: Approximation of the invariant measure n e for different values of e. These histograms are 
built from 3000 iterations of the community process Z e until time 1000. The parameters are b = 0.4, 
d = 0, c = 0.005, B = 0.02, r = 2, D = 0 and C = 0.04. 


set corresponds to the support of the stationary distribution n. 

We now compare these measures 7r £ with the measure 7f. With the parameters of the 
simulations, the averaged invariant measure 7f = l-i n b( n .h* ) satisfies 

b n ~ 1 

Mn = 7 Wfil, V n > 2 

( c) n 1 n(n!) 

OO 

'y ] Hn = 1 ) 

n= 1 


where c = c + rB 2 jC is the apparent competition. 
In the simulation, we approximate [i\ by 


1 

~ Ef,i(VJ)‘-‘KP! 


~ 2,69 • 10“ 5 . 


In Figure [9] we consider the marginal distribution of prey and predators given by the previous 
simulations (e = 1 in the left column, e = 0.1 the middle column and e = 0.00001 in the right 
column). These distributions are projections of the histograms in Figure[8j We compare them 
with the projections of the averaged distribution 7f represented with a black line joining the 
points {n,n n ). For these histograms, we chose subdivisions centered in the integers for the 
prey population and in the (h* j)>i for the predator populations. 

We observe a rapid convergence of the marginal distributions of prey and of predators, toward 
the marginal averaged distributions. 

In these simulations, the sequence ( 7 r e ) £ seems to converge in law to n £ but we have yet 
no mathematical proof of this result. 
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£ = 1 


£ = 0.1 


e = 0.00001 







Figure 9: Marginal invariant distributions of the prey population in (a 
population (d)||(e)f(f) for different values of 


)|(c) and of the predator 
EB e = 0.1 in |(b)|(e)1 and e = 0.00001 in 


e (e = 1 m 

(e)|[(f) ). These histograms are built from 3000 iterations of the community process Z s until time 1000. 
The parameters are b = 0.4, d = 0, c = 0.005, B = 0.02, r = 2, D = 0 and C = 0.04. 


5 Discussion 

We introduced new models for prey-predator communities in which the predator dynamics 
is faster than the prey one. These stochastic models derive from a microscopic model of the 
community under the successive scalings of large predator population size and small predator 
mass. In both cases, we assume that the prey population size remains finite. This assump¬ 
tion corresponds to a biological reality as in the case of forest or experimental settings: for 
example in |25j the authors study trees-insects communities in which the number of trees is 
of the order of a hundred. These models can be easily simulated and their dynamics is sim¬ 
ple therefore they represent an alternative to slow-fast Lotka-Volterra dynamical systems m- 

From an ecological point of view, we proved that our prey-predator process admits a long 
time distribution. Since the convergence to the invariant measure is exponentially fast, this 
invariant measure can be rapidly simulated using trajectories of the process. This distribution 
gives insight on the composition of prey-predator systems and can be used to compute the 
composition of trees-insects communities when the counting of insects may be difficult to do 
experimentally. Moreover, in the situation where the mass ratio between prey individuals and 
predators is small, the limiting distribution admits an explicit expression, which is valuable 
for applications. 
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In this article, our approach was to consider first a large predator population size and 
then the small size of predators. The interest of these successive scalings was to introduce 
two different models corresponding to two different biological configurations. It would also 
be possible, and biologically relevant, to consider both scalings simultaneously by considering 
accelerated birth and death events as in fl3 j (Section 4.2). In this case other processes would 
arise as the coupling of a birth death process and a diffusion. The study of such processes 
will be considered in future works. 

Knowing the long time behavior the prey-predator community is crucial to consider the 
evolutionary dynamics of the community. Indeed if we are interested in the phenotypic evolu¬ 
tion of phenotype traits of prey and/or predators (as the development of specific tree defenses 
or of specialist insects strategies [32]), one might be interested to consider the arrival of rare 
mutations as in adaptive dynamics settings (e.g. [3B, [23]). In this setting the stationary 
distribution represents the state of the resident community as a mutant arises. Therefore 
its knowledge is important to study the possible invasion of the mutant phenotype and to 
characterize the favourable strategies for prey or predators m- 

Finally, even if our motivation was to describe the dynamics of ecological communities such 
as trees-insects communities, these models or methods could also be applied to very different 
questions such as epidemiology. For example, the study of the propagation of insect trans¬ 
mitted diseases such as Malaria impose to consider the interaction between insect and human 
populations [3J. In this case the dynamics of insects is much faster than the human dynamics. 
Therefore piecewise deterministic processes could be introduce to model the dynamics of the 
disease and to study its outbreaks. 
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